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METHOD FOR CONVERTING CT DATA TO 
LINEAR ATTENUATION COEFFICIENT MAP DATA 

BACKGROUND OF THE INVENTION 
5 Field of the Invention 

[0001] This invention generally pertains to systems and methods for 
analyzing computer tomography data. More particularly, the present invention 
is directed to a method of converting computer tomography (CT) data to linear 
10 attenuation coefficient data for use in nuclear medicine, such as to 
compensate for attenuation in SPECT and PET imaging. 
Description of the Related Art 

[0002] Medical imaging falls into two distinct modalities or types. 

15 Transmission imaging refers to imaging such as X-ray imaging where the 
imaging source ( e.g, , X-ray) is external to the subject and is transmitted 
through the subject to a detector. Emission imaging refers to imaging where 
the imaging source ( e.g. , gamma-emitting radiopharmaceutical) is internal to 
the subject (as a result of injection or ingestion) and is emitted from the 

20 subject to a detector. Attenuation of source radiation occurs when the source 
radiation passes through the subject tissue, as a result of the subject 
absorbing or scattering some of the radiation photons. In general it is a 
simple matter to determine the attenuation of a discrete transmission source, 
since the amount of the external source being transmitted through the subject 

25 is known, and can be compared with the amount of radiation exiting the 

subject. However, measurement of attenuation in emission imaging is more 
difficult, because the accurate amount of emission source radiation being 
generated in the subject that results in a quantity of radiation being detected 
outside the subject cannot be measured directly. 

30 [0003] Appropriate corrections for scatter and attenuation correction 
are prerequisites for quantitative nuclear medicine. X-ray CT image volumes 
can be used to derive Linear Attenuation Coefficient (LAC) maps ("mu- 
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maps"), suitable for compensating for attenuation in single-photon-emission- 
computer-tomography (SPECT) and positron-emission-tomography (PET). 
[0004] In general, a transmission scan is performed at an energy level 
other than the energy of the emission scan. Thus, the resulting attenuation 
5 map needs to be scaled to the actual emission energy of the scan, before it 
can be used to correct for attenuation in the emission reconstruction process. 
For source-based derived mu-maps, the conversion is simple because the 
discrete transmission and emission energies are known. For x-ray CT 
however, the transmission spectrum is continuous (and not discrete as it is 
10 the case for source-based methods), and, more importantly, depends upon 
the particular CT scanner and the attenuating body. 

[0005] Figure 1 shows that attenuation coefficients for different types of 

tissue depend on the energy of the photons, and can be grouped in 
essentially four groups, depending on their atomic number, Z: Air, soft tissue, 
15 bone, and iron, with iron representing a class of "Very High-Z" implants, such 
as surgical screws, hip-replacements, or other possible very high-Z materials 
in the body. 

[0006] X-ray CT images are calibrated so that each voxel is measured 
in units of Hounsfield, usually defined as: 

20 HU Material = ^Material _ ^Waterj , ^.Water _ ^Air^ * ^qqq (1 ) 

. In this definition p T Matenal j s the linear attenuation coefficient of a given 
material at an "effective" transmission energy T, and p T Water is the linear 
attenuation coefficient of water at the same "effective" transmission energy T. 
The linear attenuation coefficients are "narrow beam" values, which are 
25 derived from primary photon counts only, and thus do not include any 
scattered photons (see also Fig. 16). 

[0007] Because a CT scanner emits a continuous spectrum of x-rays, 
an "effective" transmission energy is not easily obtained. The HU values are 
normalized for all scanners and protocols if the CT scanner for clinical 
30 practice has been properly set-up, so that water always corresponds to HU=0 
and air corresponds to HU= -1000 (it is noted that some definitions multiply by 
a factor of 2 10 = 1024; such cases are included within the scope of the 
invention). All clinical CT scanners have to be calibrated using a vendor 
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specific protocol to conform to this definition. However, there is no definition 
for densities greater than water. For instance, the same bone tissue may have 
different HU values when acquired with different CT scanners. The HU value 
of a bone specimen may even change depending on the surrounding amount 
5 of soft tissue and reconstruction parameters on the same CT scanner. 
Converting bone tissue accurately is important because it may contribute 
largely to attenuation of emission energy. 

[0008] Various approaches are known for converting CT values to 
linear attenuation coefficients, depending on the degree and type of 
10 approximation treating the continuous CT spectrum. These methods can be 
grouped in two major classes: Uniform scaling, and Non-uniform scaling. 
[0009] In Uniform scaling, all pixels in the transmission slices are 

multiplied by the same factor K, where K is usually computed from the linear 

_ Water 

attenuation coefficient of water at the effective transmission-energy T Ut 

__ Water 

15 and emission-energy E p E 

K = M E Water /MT Water (2). 
[0010] This approach is accurate for soft tissues (low Z), since their 

attenuation properties are similar to water. In water, or other low-Z materials, 
Compton scattering is the dominant effect for clinically used emission 

20 energies, but for high-Z materials such as bone, the photoelectric absorption 
becomes the dominant effect. As a result, a scaling factor derived using 
water deviates considerably from a bone-derived scaling value. Thus, 
equation (2) provides inaccurate linear attenuation coefficients for materials 
more dense than water and soft tissue. The equation also assumes that an 

25 "effective" transmission-energy is known. 

[0011] For a CT scan, the transmission energy can only be an 

"effective" transmission-energy, which essentially replaces the continuous 
transmission spectrum through the body with a mono-energetic "effective" 
transmission-energy. However, such an effective transmission-energy 

30 depends on the type of CT scanner, patient, and also on the CT 

reconstruction parameters. Examples of such methods are described in Kalki 
K., Brown, J.K.,et al., "Myocardial Perfusion Imaging with a Correlated X-ray 
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CT and SPECT System: An Animal Study", IEEE TNS, 43(3), 1996, pp. 2000- 
2007, 1996, which is incorporated herein by reference. 
[0012] Non-Uniform Scaling Tissue method allows regions of the CT 
image volume to be defined as part of different classes of tissues 
5 ("segmentation"). Either each class of tissue is scaled from some "effective" 
transmission energy to the emission energy, or pixel values for each tissue 
type are replaced with the appropriate attenuation coefficients at the emission 
energy. Typical choices for tissue types include soft tissue, bone, lung and air. 
However, there are numerous problems with this approach. For example, this 
10 method does not take into account existing variations in tissues densities for 
the same tissue classes, it is limited to the errors of the segmentation routine, 
it may require an educated user to segment the data accurately, and it is not 
user-friendly. 

[0013] Pixel-by-pixel conversion is an extension of the tissue typing 
15 approach. In this method, each pixel is scaled from CT units to linear 

attenuation coefficients. In principle, this method requires knowledge of the 
type of tissue for each pixel. This is often difficult because pixels may contain 
more than one tissue type ("partial volume"), or an educated user may be 
needed to identify the tissue type of each pixel (see, D.R. White, M. 
20 Fitzgerald, "Calculated attenuation and Energy Absorption Coefficients for 
ICRP Reference Man (1975) Organs and Tissues", Health Physics, 33, pp. 
73-81, 1977, which is incorporated herein by reference). Therefore, this 
method suffers from similar problems. 

[0014] Thus, there is a need for new and improved systems and 
25 methods for creating linear attenuation coefficient maps from CT image data 
that improve upon and solve the foregoing problems with the prior art. 

SUMMARY OF THE INVENTION 

[0015] In accordance with the foregoing and other objects, the present 
invention provides systems and methods for converting CT data to a linear 
30 attenuation coefficient map for use in nuclear medicine. 

[0016] According to one embodiment of the present invention, a simple 
method is provided for converting CT images obtained from an arbitrary 
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clinical CT scanner to "mu-maps". This method adapts to patient specific data 
and the varying parameters of the CT scan and reconstruction, and thus 
eliminates the need for any additional calibrations, other than the clinically 
necessary, vendor specific CT scanner calibrations. The present invention 
5 represents an optimal trade-off between accuracy and simplicity. 
[0017] The present invention provides a number of non-limiting 
advantages. For example, the present invention is capable of the following: 

converting to emission energies from 50 keV to 51 1 keV within a 
similar accuracy range; 
10 - modeling the energy dependence of tissues accurately over the whole 
clinically used emission energy range (50, 511) keV; 

estimating from the clinical CT images an effective energy; 

converting to an effective emission energy in the case of multiple 
emission peaks; and 
1 5 - adjusting the linear attenuation map to the acquisition window size in 
the emission scan. 

[0018] According to an embodiment of the present invention, a method 
is provided for converting output from a computer tomography (CT) device to 
a linear attenuation coefficient map. The method includes the steps of 

20 receiving output pixel data from a CT device for a pixel of a CT image; 

comparing a value of the pixel data to a predetermined range; if the value is 
within the predetermined range, calculating a linear attenuation coefficient 
from the pixel data using a first function; if the value is outside the 
predetermined range, calculating the linear attenuation coefficient from the 

25 pixel data using a second function; and storing said calculated coefficient in a 
memory as part of a linear attenuation coefficient map. 
[0019] According to an embodiment of the present invention, a method 
is provided for converting output from a computer tomography (CT) device to 
a linear attenuation coefficient map. The method includes the steps of: 

30 receiving output pixel data from a CT device for each pixel of a CT image and 
for each pixel: comparing a value of the pixel data to a predetermined range; if 
the value is within the predetermined range, generating a first function and 
calculating a linear attenuation coefficient from the pixel data using the first 
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function; and if the value is outside the predetermined range, generating a 
second function and calculating a linear attenuation coefficient from the pixel 
data using the second function; and generating a linear attenuation coefficient 
map based upon each linear attenuation coefficient calculated for each pixel. 
5 [0020] According to an embodiment of the present invention, a system 
is provided for generating a linear attenuation coefficient map. The system 
includes an input, a processing unit and an output. The input is for receiving 
pixel data. The processing unit is configured to compare a value of the pixel 
data to a predetermined range, to select a function based on whether the 
10 value is within the predetermined range or outside the predetermined range, 
and to calculate a linear attenuation coefficient from the pixel data using the 
selected function. The output is for outputting the calculated linear 
attenuation coefficient. 

[0021] According to an embodiment of the present invention, a 

15 computer program product is provided which resides on a computer readable 
medium, for interactively constructing, editing, rendering and manipulating 
pixel data for a CT image. The computer program includes a computer 
executable instructions for causing the computer to perform the following: 
receiving output pixel data from a CT device for a pixel of a CT image; 

20 comparing a value of the pixel data to a predetermined range; 

if the value is within the predetermined range, calculating a linear attenuation 
coefficient from the pixel data using a first function; if the value is outside the 
predetermined range, calculating the linear attenuation coefficient from the 
pixel data using a second function; and storing the calculated coefficient in a 

25 memory as part of a linear attenuation coefficient map. 

[0022] With these and other objects, advantages and features of the 
invention that may become hereinafter apparent, the nature of the invention 
may be more clearly understood by reference to the following detailed 
description of the invention, the appended claims, and the drawings attached 

30 hereto. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0023] The invention will be described in detail with reference to the 

following drawings, in which like features are represented by common 
reference numbers and in which: 

[0024] Fig. 1 shows linear attenuation coefficients for different types of 

tissue and their dependencies on the energy of the photons; 
[0025] Fig. 2 is a typical interpolated graph demonstrating the 

relationship between CTs Hounsfield Units (HU) and linear attenuation 
coefficients; 

[0026] Fig. 3 is a typical conversion graph from HU to linear attenuation 

coefficients at 100 keV transmission energy and 150 keV emission energy; 
[0027] Fig. 4 is a graph showing the most prominent interaction for 

biological materials and relevant energies; 

[0028] Fig. 5 is a conversion graph from HU to linear attenuation 

coefficients approximated as piecewise linear functions; 

[0029] Fig. 6 is a flow diagram of a method for converting from HU to 

linear attenuation coefficients according to an embodiment of the present 

invention; 

[0030] Fig. 7 is a graph of soft tissue normalized attenuation 

coefficients; 

[0031] Figs. 8a-8b are a slice of CT of a torso region and a 

corresponding histogram; 

[0032] Fig. 9 is a histogram of the CT volume of Fig. 8b at a resolution 

of 500 bins; 

[0033] Figs. 10a-10b are a slice of CT of a pelvis region and a 

corresponding histogram; 

[0034] Figs. 1 1a-1 1 b are a slice of CT of a skull region and a 

corresponding histogram; 

[0035] Fig. 12 is a flow diagram of a method for computing the effective 

transmission energy, according to an embodiment of the present invention; 
[0036] Fig. 1 3 is a chart of the optimal default effective transmission 

energy, according to an embodiment of the present invention; 
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[0037] Fig. 14 is a chart showing error estimation of the effective 
energy, according to an embodiment of the present invention; 
[0038] Fig. 1 5 is a diagram of a collimator showing "narrow beam" 
attenuation of initial flux; and 
5 [0039] Fig. 16 is a diagram of a collimator showing "broad beam" 

attenuation of initial flux. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
[0040] A typical interpolated graph demonstrating the relationship 

between CT's Hounsfield Units (HU) and linear attenuation coefficients is 

10 shown in Fig. 2. The conversion from HU to linear attenuation coefficients is 
dependent upon emission-energy and transmission-energy. For instance, a CT 
image acquired at two different effective transmission energies ("fx"), will need 
to be scaled differently to the same emission energy, as can be seen by 
comparing Fig. 2 and Fig. 3. This difference in scaling occurs because of the 

15 different physical effects that contribute for varying energies at various 

amounts to the attenuation. For energies below 511 keV, the attenuation can 
be considered as a combination of both Compton scattering and photoelectric 
absorption. But the higher the energy, the lower the contribution of the 
photoelectric absorption is to the attenuation, and at 51 1 keV, it is essentially 

20 negligible for all biological materials (including bone). In the CT energy range 
(50-150 keV), however, the photoelectric absorption is not negligible (see Fig. 
4 and also see, Sorensen J.A., Phelps M.E., "The Physics of Nuclear 
Medicine", Grune & Stratton, 1980, which is incorporated herein by reference). 
[0041] The conversion mechanism used by the present invention is 

25 based on the fact that some types of media are always mapped onto certain 
HU values. 

Air - -1000 HU independent of effective transmission 

energy 

Water - 0 HU independent of effective transmission 

30 energy 

Cortical Bone - - 1000-2000 HU dependent on effective transmission 
energy 
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[0042] The scaling between the energy dependent linear attenuation 

coefficient and the HU, which are essentially water normalized differences 
between the attenuating material and water, can be approximated as piecewise 
linear functions (see Fig. 5), for the transition between low Z (soft tissue) and 
5 higher Z (bone) materials. For water, the linearity is exactly based on the 
definition of the HU. 

[0043] For the purposes of the present invention, and without loss of 

generality, two regions shall be considered: air-to-soft tissue(water) (HU<=0); 
and soft tissue(water)-to-bone(HU>0). 

10 [0044] Fig. 6 is a flow diagram of a method of creating a mu-map 

according to the present invention. It will be understood by those skilled in the 
art that the following processes may be programmed onto central processor or 
computer, which may be coupled with a CT scanning device, or which may 
receive CT data in any fashion. Accordingly, a program written to perform the 

15 following may be stored conventionally in memory. Accordingly, the present 
invention is not limited to any particular system configuration. 
[0045] Beginning with step 6-1 , CT pixel values are converted to HU 

values. CT devices may output pixel data a number of ways. Typically, a 
device outputs data in a proprietary fashion that must be converted to HU. The 

20 conversion is usually provided by the vendor. For example, CT devices may 
output DICOM ("Digital Imaging and communication in Medicine" standard) 
compliant data containing a header that defines the conversion from 
proprietary pixel data to HU. Obviously, if a CT device outputs pixel values as 
HU values, then the conversion is unity. 

25 [0046] At step 6-2, HU (bone) and the effective transmission energy are 

estimated. In one embodiment of the present invention, a simple "numerical 
observer" replaces the educated human observer. However, it is to be noted 
that the present invention also applies to more sophisticated numerical 
observers, using ejj. Artificial intelligence, data base scoring, statistical 

30 methods, etc. One implementation of such a "numerical observer" uses the 
histogram of the CT image volume to find the most common bone HU value at 
the given scan parameters in the CT image volume. This value is assigned to 
be HU T Bone , where the sub-script "T" indicates that this value is at some 
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"effective" CT transmission energy. This "effective" CT transmission energy 
includes all possible effects that could influence the bone-HU, including effects 
from, e.g., the CT reconstruction filter. 

[0047] A flow diagram of the "numerical observer" is shown in Fig. 12. 

5 At step 12-1 , the histogram is generated from the CT data for the entire CT 
image, by conventional means. The histogram is binned by HU values and 
number of pixels per bin. The histogram of a CT image volume of a patient is 
normally dominated by soft tissue. Bone values span a wide range and the 
abundance thereof is, in general, orders of magnitude less than the abundance 

10 of soft tissue values. Hard cortical bone is rather rare and can be found mostly 
as a shell around the bone structure. Furthermore, bone density decreases 
with age and is smaller for adolescents than for adults. Thus, it is a real 
challenge to find the bone value most suitable to set the conversion function 
from HU to linear attenuation coefficients. 

15 [0048] It is assumed that the bone-histogram of a CT image volume of 

humans can fall in two distinct classes, if the histogram is properly binned: (a) 
one dominant bone peak (e.g. Head) in the histogram, or (b) no dominant bone 
peak (e.g. Pelvis). 

[0049] The HU values of bone are clearly larger than zero (0) because 

20 the linear attenuation coefficient of bone is, in most cases and for the relevant 
energies, roughly twice as large as that of water and soft tissue (see Fig. 7; the 
linear attenuation of bone is roughly 1 .75-2 times larger than of water and soft 
tissue). It is therefore possible to search for the bone peak above a threshold 
HU-value. The binning parameters of the histogram need to be adjusted to 
25 obtain a slowly varying function to search for HU T Bone , yet it needs to be 

sensitive enough to distinguish between the main sections of the human body 
with varying amount of bone tissue. 

[0050] The optimal CT image histogram is preferably binned in N bins 

in the range (HU B0 ne,Min, HU Bo ne,Max) where the maximum value of the histogram 
30 is selected as estimate for HUt 8006 . The range (HU Bon e,Min, HUeone.Max) 
corresponds to a CT effective energy range of and is chosen to include all 
clinical X-ray CT scanners. Other physical effects, such as reconstruction filter 
choice, are essentially included into this value. 
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[0051] At step 12-2, the location of the local maximum value (dominant 

bone peak) is determined. If there is a dominant bone peak, then the location 
of the dominant bone peak is used as the estimated HUT bone . If there is no 
local maximum, then at step 12-3 the Laplacian (i.e., 2 nd derivative) of the 
5 histogram is generated. Then, at step 12-4, the location of the local maximum 
value of the Laplacian is determined. If there is a local maximum, it is used as 
the estimated HU T bone . 

[0052] As examples, Fig. 8b and 8a, respectively show a slice of the 

torso region (Kidney) and the histogram, binned with parameters. The lower 

10 curves show first and second derivatives. In this case a peak is clearly present 
representing a range of HU value, corresponding to the most prevalent bone 
structure. The kidney stone is rather small and localized and is represented as 
the second peak at higher HU units, which is not used to derive the effective 
energy. Fig. 9 shows the histogram with a rather fine binning, showing the 

15 values for bone vary strongly. These figures shows how 'noisy" the data is. 
The lower curves show first and second derivatives. It would be very difficult to 
find the most relevant bone value from this histogram. Figs. 10a-10b and 11a- 
11b show slice and histograms of other regions of the body, representing the 
most common bone value distributions, which can be analyzed with the 

20 "numerical observer" of the present invention. In Figs. 10a-10b, The lower 

curve shows the second derivative. The location of the maximum of the second 
derivative is used to estimate the highest reasonable bone value. No distinct 
bone peak. In Figs. 11a-11b, note the distinct, but wide, bone peak from the 
skull region. The lower curves show the first and second derivatives. 

25 [0053] Once the most likely value of bone tissue HU T Bone in the CT 

image volume is estimated, this value can be used to estimate the "effective" 
CT transmission energy using equations (1) and (2). First, equations (1) and 
(2) are used to find (j E Bone corresponding to the estimated HU T Bone as follows: 
HU T Bone = (MT Bone -Mi Water ) / (Mr Water -Mr Ajr ) * 1000; (3) 

30 ME Bone / MT Bone = ( a*. * + c^e) / ( a^ * T^ + Cbone ) = K Bone (E, T). 

(4) 

[0054] Here (3) is the Hounsfield expression for bone tissue, a^ , b bone , 

and c c bone are bone coefficients from a fit, and E and T are user-defined 

11 
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emission and estimated transmission energies. Combining (3) and (4), |j E Bone is 
expressed through HU T Bone as: 

Me Bone = K Bone ( E> T) * (( HUj **» / 1000 )*( ^Water _ ^Air ) + ^Water y (5) 

[0055] The effective transmission energy T is solved from (jE Bone with a, 

5 b, c being known for bone and water. 

[0056] Even if no appropriate bone value is found, the CT effective 

energy needs to be determined. Accordingly, if there is no local maximum, then 
a default effective transmission energy is selected at 12-6. An error analysis 
allows the computation of the optimal default energy, which minimizes the error 
10 of the conversion (see Fig. 13). 

[0057] Once the effective transmission energy is determined at step 6- 

2, conversion functions F and G are created at step 6-3. The energy 
dependence of the linear attenuation coefficient can be fitted by a double- 
power law, with coefficients a, b, c for energies below and above a transition 

15 energy E T h*. 

p(E) = a<E b< + c< for E < Eth; 
p(E) = a>E b > + c> for E > Eth. (2) 
Above the transition energy, Compton scattering is much more prevalent than 
photoelectric absorption for low-Z materials. The coefficients a, b, c are 
20 "weighted" X 2 fit to data (e.g., data published by J.H. Hubbell (see references 
below) and can be obtained from the National Institute of Standards and 
Technology, e.g., see www.nist.gov ). 

[0058] The X 2 fit is done for below 150 keV, and above 150 KeV, 

producing two sets of a, b, c values. The X 2 fit is "weighted" but not in the usual 

25 manner. Weights are typically derived from the variance of the error of a 

particular data point. In this case, weights have the effect that data points with 
larger measurement error don't "contribute" as much to the fit as data with a 
small error. This is used to "nail" the fit at the transition energy better, which in 
turn gives a fit which is a) continuous in the first derivative for all E in the 

30 energy range, and b) has very small deviation from all data points, but in 

particular at the clinically important energy region around 150keV. The ratio of 
the weight at 150keV with respect to all other weights at the other energy 
values is 1/10. 

12 
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[0059] The main contribution to the attenuation coefficient comes from 

the varying cross sections of the photoelectric and Compton effects, which in 
turn also scale with Z x and E y (where x, and y depend on the effect) changes 
around the chosen transition energy E t h = 150 keV. This transition energy is 
5 chosen in part, because in nuclear medicine imaging, 150 keV is close to 
140keV - the TC99 emission energy (Tc99 is very often used in SPECT). 
[0060] For iron (Fe), the transition energy selected to be 250 keV, 

which is also fitted in a similar fashion. Although iron is a material that does not 
normally occur in the body, it is still found in many people (e.g., surgical nails, 
10 implants, bullets, wires, etc.), and thus, there is a need to properly compute the 
linear attenuation coefficient. 

[0061] The errors of the fit for the mentioned soft tissues and water < 

0.2% and for bone <0.55% for all energies. The function a*E b +c is very easy to 
use, in the subsequent calculations, and is well suited for the purpose. 
15 [0062] At step 6-4, it is determined whether the pixel value is between - 

1000 (HU for air) and 0 (HU for water). If so, processing proceeds to step 6-5. 
If not, processing proceeds to step 6-6. 

[0063] At step 6-5, air-soft tissue region is converted using a linear 

function F, calculated from the data points of air and water, for which the 
20 energy dependence is obtained from the above fit to the double-power law (2). 
Because air and water values on the HU scale are independent of the 
transmission energy, this conversion is accurate. 

[0064] Soft tissue-Bone approximation (HU>0) is more difficult. The HU 

value of bone is scanner, scan, reconstruction, and patient dependent and 
25 depends on the effective transmission energy, which is generally unknown. A 
human observer could identify the bone tissue in a CT volume and assign 
linear attenuation coefficients to this type of tissue, which is indexed by the HU 
values. 

[0065] At step 6-6, a conversion for all other values can then be 

30 calculated using the conversion function G for the Soft tissue-Bone region. 
Function G is approximated by a preferably linear fit to the data from soft-tissue 
and the identified bone tissue, with the constraint that the conversion function 
is continuous over the whole HU-range, i.e. F(HU=0)=G(HU=0). 

13 
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[0066] Thus, values for |j are generated from the pixel values of any x- 

ray CT imaging device. The generated p-map may be used in nuclear 
medicine imagining, such as SPECT and PET. 
[0067] Estimated Errors 

5 [0068] Because both F and G functions are pixel value dependant, it 

makes this approach more accurate than the Tissue Type approach. On the 
other hand, the fact that two different conversion regions are considered 
separately makes it more accurate than Uniform Scaling approach. The main 
limitation of this algorithm, besides small errors caused by linear interpolation 

10 of F and G, is due to the uncertainty in CT effective energy for the Water-Bone 
region. Fig. 14 shows the whole range of possible mis-estimation of the 
effective transmission energy and the resulting error in the LAC. In general, the 
effective transmission energy can be estimated within ±10%, and thus can 
expect the error in the LAC to be roughly ±5%. This can, of course, be 

15 improved by calibration measurements. 

[0069] Multiple Energy Groups for Emission Data 

[0070] In nuclear medicine application SPECT acquisition may be 

performed with multiple acquisition windows summed into one projection 
image. The conversion problem can be generalized to a conversion from the 

20 effective transmission energy to the effective emission energy, wherein a 
discrete spectrum can be replaced by one "effective" emission energy value. 
[0071] The method to calculate the appropriate effective emission 

energy can either use a branching-ratio averaged energy or more accurately a 
weighted X 2 fit, which can include n windows. Using the branching-ratio 

25 averaged energy method: if the emission energies are Ej, with i = 1 n, and 

branching ratios (or mixing ratios) br it then the "effective" emission energy can 
be computed using the following equation: 

Eeff.emission = Ei*bri +...+ En*br n , where 
bn +...+ br n =1. 

30 [0072] This formula is only an approximation, but the errors are usually 

less than 1%. A more accurate method is obtained by X 2 -fit having fit 
parameters that are the effective narrow beam linear attenuation coefficient 
and un-attenuated counts. The X 2 -fit results are more accurate and preferred. 

14 
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[0073] Comp nsation for Acquisition Window Width 

[0074] In practice, SPECT projection data always contains object 

scatter, due to necessary large size of the acquisition window, and thus, 
compensating for attenuation in a reconstructing with those "narrow beam" 
5 attenuation values, may result in an over-correction. A common method to 
correct for this problem is to reduce the attenuation values by some factor. This 
factor depends on the energy of the nuclear medicine image and the 
acquisition window width (the wider acquisition window the more scatter is 
detected), which is often neglected. These correction factors are computed 
10 from measurements, and a X 2 -fit can be used to adjust for all possible energies 
and window widths. 

[0075] Attenuation is the effect of detecting less photon flux due to the 

presence of an attenuating medium. If no medium is present the initial flux is 
denoted in the following as lo and in the presence of a material it is denoted as 

15 l(x, p), where x is the length of the traversed distance of the photons in the 
material with linear attenuation coefficient p. The value l(x, p) and lo are related 
by Beer's law: / (x,fj) = \oe~ yx (see Fig. 15). The linear attenuation coefficient p 
in this equation is assumed not to contain any scatter, and is referred to as the 
"narrow beam" linear attenuation coefficient. 

20 [0076] The difference between "narrow beam" and "broad beam" 

coefficients is that the first does not include any object scatter, and is the ideal 
case, the latter does. The name refers to the experimental set-up, where the 
beam is highly collimated before entering the detector, and thus is a narrow 
beam, or where it is not collimated and thus is a broad beam (see Fig. 16). 

25 Object scatter is allowed to enter the detector, which could be accounted for in 
an object-dependent build-up factor, not discussed here. 
[0077] The linear attenuation coefficient p is a material property 

quantifying ability to attenuate radiation. A material of thickness 1 cm and linear 
attenuation coefficient p 1 1 reduces the initial flux from l 0 to l 0 /e M * 1cm . The unit 

30 of the linear attenuation coefficient is 1/cm. The linear attenuation coefficient is 
energy dependent, where higher energy radiation is more penetrating than 
lower energy radiation, thus the linear attenuation for a given material is larger 
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(more attenuating) at lower energies than at higher energies. Example: 
Pwater(140 keV) = 0.15847 1/cm < Mwater(70 keV) = 0.19275 1/cm. 
[0078] Thus, the invention has been described in connection with what 

are presently considered to be the most practical and preferred embodiments. 
5 However, it is to be understood that the invention is not to be limited to the 
disclosed embodiments, but, on the contrary, is intended to cover various 
modifications and equivalent arrangements included within the spirit and scope 
of the appended claims. 

[0079] For example, it should be noted that the present invention may 

10 be extended to include numerous additional functionality. For example, 
contrast agents may be handles in a conversion from CT data to linear 
attenuation coefficients. Other numerical observers may be used to improve 
the estimation of the effective energy using e.g., CT parameters from DICOM 
header, Template matching of the CT image volume to a reference CT image 

15 of the scanner, Extending the numerical observer to other modalities, such as 
MRI. Compensation for the acquisition windows width may use the computed 
narrow beam mu-map to compute a "scatter projection estimate" for the actual 
acquisition window settings. This scatter projection estimate has the advantage 
that it is patient dependent and non-local. Mu-maps may be used to 

20 classify/segment tissue/voxel types accurately. The results of segmentation 
may be used for automated analysis of co-registered images from other 
modalities, e.g. MRI, Angiography, or Ultrasound. The results of segmentation 
may be used for automatic registration with atlas CT images, improved mutual 
information registration. Quantitative analysis may be performed on the data 

25 and displayed. 

[0080] For further information, the following references, each of which 

is incorporated by reference, provide background information and/or details the 
will assist the reader in understanding the present invention: 

Bai C, Shao L, Da Silva A.J. et al., A Generalized Model for the Conversion 
30 from CT Numbers to Linear Attenuation Coefficients", IEEE MIC 2002; 

Blankespoor S.C., Wu X., et al., "Attenuation Correction of SPECT Using X- 

ray CT on an Emission Transmission CT System: Myocardial Perfusion 

Assessment", IEEE TNS, 43(4), 1996, pp. 2263-2274; 
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Kinahan P.E., Townsend D.W., et al., "Attenuation Correction for a combined 
3D PET/CT scanner", Med. Phys. 25(10), 1998, pp. 2046-2053; 
Kalki K., Brown J.K., et al., "Myocardial Perfusion Imaging with a Correlated 
X-ray CT and SPECT System: An Animal Study:, IEEE TNS, 43(3), 1996, pp. 
5 2000-2007; 

E.B. Saloman, J.H. Hubbell, and J.H. Scofield, "X-ray Attenuation Cross 
Sections for Energies 100 eV to 100 keV and Elements Z = 1 to Z - 92", 
Atomic Data and Nucl. Data Tables 28, 1 (1988); 
Hubbell, J.H. (1982), Photon Mass Attenuation and Energy-Absorption 

10 Coefficients from 1 keV to 20 MeV, Int. J. Appl. Radiat. Isot. 33, 1269-1290; 
Hubbell, J.H. (1994), Bibliography of Photon Total Cross Section (Attenuation 
Coefficient) Measurements 10 eV to 13.5 GeV, 1907-1993, NISTIR 5437; and 
Hubbell, J.H. Gerstenberg, H.M. and Saloman, E.B., Bibliography of Photon 
Total Cross Section (Attenuation Coefficient) Measurements 10 eV to 13.5 

15 GeV, National Bureau of Standards Internal Report NBSIR 86-3461 (1986), 
which is incorporated herein by reference. 
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